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Abstract 

We discuss 77 — )• tttt reactions, starting from the two-pion threshold up to about ^/s^ = 6 GeV. 
Several reaction mechanisms are identified. We include the dipion continuum due to pion exchange 
(for 7r+7r~) and exchange (for vr'^vr'') as well as the pronounced dipion (j(600), /o(980) and 
/2(1270) s-channel resonances. We discuss also a possible contribution of less pronounced scalar 
resonances /o(1500) and /o(1710) being glueball candidates, as well as the tensor the resonances 
72(1525), /2(1565) and /2(1950). While the contribution of /2(1525) is rather small, the contribu- 
tion of /2(1565) changes the spectrum around the resonance position. We find that the inclusion 
of the spin- four /4(2050) resonant state improves the situation at ^/s^ « 2 GeV. We estimate the 
relevant diphoton partial decay width to be T f^^^^ ^ 0.7 keV. At higher energies, Brodsky-Lepage 
and hand-bag mechanisms are included in addition. We nicely describe the world data for 77 — )• tttt 
for the first time both for the total cross section and angular distributions for 77 — )• vt^tt^ and 
77 — )• vr'^vr'^ reactions simultaneously at all experimentally available energies. 

The cross section for the production of two-pions in ultraperipheral ultrarelativistic heavy ion 
collisions is calculated in the impact parameter space Equivalent Photon Approximation (EPA). 
We obtain the nuclear cross section of 46.7 mb for ■k~^-k~ and 8.7 mb for vr'^vr'^ at the LHC en- 
ergy y/SMN = 3.5 TeV and minimal cut on pion transverse momentum pt > 0.2 GeV. Differential 
distributions in impact parameter, dipion invariant mass, single pion and dipion rapidity, pion 
transverse momentum as well as pion pseudorapidity are shown. The 77 — t- tt+tt^ subprocess con- 
stitutes a background to the exclusive AA — )• Ap^{-^ tt'^'k~)A process, initiated by photon-pomeron 
or pomeron-photon subprocesses. We find that only a part of the dipion invariant mass spectrum 
associated with 77-collisions can be potentially visible as the cross section for the AA — )• Ap^A 
reaction is very large. 

PACS numbers: 25.75.-q, (Heavy-ion nuclear reactions relativistic) 
25.75.Dw, (Particle production (relativistic collisions)) 
24.30.-v, (Nuclear reactions: resonance reactions) 
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I. INTRODUCTION 



Ultrarelativistic ultraperipheral production of mesons and elementary particles is a special 
category of nuclear reactions The ultrarelativistic ions provide large fluxes of quasi-real 
photons which can collide leading to different final states. Both total photon-photon cross 
sections as well as cross section for particular simple final states are interesting. The nuclear 
cross section is usually calculated in the Equivalent Photon Approximation in momentum 
space ji, S] , or in impact parameter space jl, 0] • In the latter case the flux of photons depends 
on the transverse distance from the heavy ion trajectory. Impact parameter space is very 
convenient to exclude cases when both heavy ions collide, i.e., when they do not survive the 
high-energy collision. In the past we have performed also full momentum space calculation 
for /i"*"//" production Q. In the momentum space calculation (EPA or full calculation) the 
effects of nucleus-nucleus collisions and their associated break-up are neglected. This effect 
is rather small for light particle production such as e~^e~ often studied in the literature. 

Recently we have studied several processes initiated by the photon-photon collisions such 
as: p°p° (sf, 0], QQ 0], DD 0] and for high mass dipions [8]. We have shown there 

that the inclusion of realistic charge form factors, being Fourier transforms of realistic charge 
distributions, is crucial for estimating reliably the nuclear cross sections. 

Here, we shall apply the previously used method for the exclusive production of tt+tt" and 
7r°7r°. The present analysis is an extension of the analysis of high- mass dipion production 
[sl where we concentrated exclusively on perturbative QCD effects. 

The elementary reactions 77 — > tt^tt^ and 77 — t- tt+tt^ are interesting by themselves, 
since understanding the mechanism of the 77 — )■ tttt reaction at low or intermediate energies 
{J's^ < 1.5 GeVHs very important for applications of chiral perturbation theory and pion- 



pion interaction_ 9Nl4| . At higher photon-photon energies the Brodsky-Lepage [8|, llSMlS 
and hand-bag |19|] mechanisms have been discussed in the literature. 

In addition, we study the contribution of glueballs or glueball candidates [20|. Only esti- 
mates of the total cross section were presented in the literature 2l| recently. In general, it 
is not possible to observe glueballs, even if the corresponding cross section is sizeable. One 
should concentrate on a concrete final state. The vrvr channels are a good example. Estima- 
tion of the background and its interplay with resonant signal should be better understood. 
Here we show how the glueball candidates can be seen in reality when analyzing e.g. dipion 
spectra. 

In the present paper we construct a hybrid, multi-component model for the 77 tttt 
reactions which describes all experimentally available world data. We shall extend the range 
of the photon-photon energies to the region where the transition from nonperturbative QCD 
to pQCD may be expected. It is very inaportant to understand the onset of the pQCD effects. 
For example the hand-bag mechanism [l9| is usually fitted to experimental data for energies 
bigger than 2.5 GeV. 

The 77 — ?► TTTT processes are fairly complicated. Different mechanisms may contribute 
in general. The present analysis is partially an update of an old analysis in Ref. M. In 
the present paper we shall try to understand both 77 — )■ tt+tt" as well as 77 — )■ tt^tt^ 
processes simultaneously, from the kinematical threshold {W = 2m.^) up to the maximal 
experimentally available energy W.y^ ~ 6 GeV. We will include both soft and hard continuum 
processes as well as s-channel resonances. Here we shall include many more resonances 
compared to Ref. 1611. I n addition, we shall try to describe the new very precise data of the 



Belle Collaboration [22], including the angular distributions of pions. 
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Our paper concentrates on realistic predictions of the cross sections for dipion produc- 
tion in ultraperipheral ultrarelativistic heavy ion colhsions. We shall present first realistic 
predictions for PbPb — )■ Ph{Ti'K)Ph reaction at the LHC Both total cross sections and dif- 
ferential distributions will be shown and discussed for the first time in the literature. The 
photon-photon induced dipion production in nuclear collision, interesting by itself, consti- 
tutes a background to another type of nucleus-nucleus reactions induced by phot on-p omeron 

and its 



(pomeron-photon) exchanges leading to a coherent production of meson |23|-|25 
radial excitations. The interplay of the both processes was not discussed so far in the 
literature. 



II. MODELLING 77 tttt REACTIONS 
A. The 77 vr+vr" continuum 

The Born amplitude for the 77 — )■ tt+tt^ process is shown in FigJT] The helicity dependent 
^AAAA>^ - > ^AAAA>^^ V, 

^AAAA>'- - > - - ^AAAA/ ^ ^ 

TT TT TT 

FIG. 1: The Born term matrix elements for non resonant pion pair production. 



amplitude for point-like pions is a sum of the three terms (see e.g. |26l| ) : 
contact amplitude: 

Ml[\^ (gi,g2,P.+ ,P.-) = e'J22g'^''e,i\^)e,i\2) , (2.1) 
t-channel pion-exchange amplitude: 

Ml[\^ (gi, g2,P.+ ,P.-) = E C^Pn- " Qi) ^Mi^ ^1) {^pI+ - ^2) ^^2. ^2)— ^ (2.2) 
and u-channel pion-exchange amplitude: 



u — ml 



The full pion-exchange amplitude is a sum of contact + t + u amplitudes 

MU, = Ml[\^ + Ml'^,^ + A^I;",, . (2.4) 
Such an amplitude can be formally written as: 

Ml,.^ = £^(gi, \i)e,{q2, \2)M^'' • (2.5) 
It is straightforward to show that 

ga^-M"'^" = , (2.6) 
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which is consistent with gauge invariance. The QED Born amphtude for the 77 — t- vr+vr" 
reaction with point-hke particles has been known for a long time. An interesting problem 
is to construct the QED amplitude for real, finite-size, pions. We use an idea proposed by 
Poppe jl?! , and correct the QED amplitude by an overall t and u dependent form factor: 

Ml^,^ (t, u,s) = n {t, u, s) Mlf^"" (t, u, s) . (2.7) 

It is natural that finite size corrections damp the QED amplitude for both t and u large. At 
sufficiently high energies the following simplification fulfils this requirement: 

■ (-) 

In practice one can take F{x) = exp(^^a;); = (4-6) GeV~^. F{x) denotes the standard 
vertex function which provides the convenient normalization F {0) = 1 and additionally 
F (t) — )■ when t — )• —00. In the limit of large s: 

n{t,u,s) ^ F^t) , 

n{t,u,s) ^ F^u) . (2.9) 
This s limit generates standard vertex form factors and additionally at z = cos 6 ~ 0: 



B. s-channel 77 — t- resonances 



In Figj2] we show a list of several possible dipion resonances that could, in principle, 
contribute to the 77 tctt processes. 




/o(600) 
/o(980) 

/2(1270) I 
/o(1500)^ ,.r 

/2(1565) ; 
/o(1710) 
/2(1950) 
/4(2050) 



FIG. 2: 77 — )• resonances — )• vr"'"' vr 



+/0-rr-/0 



In Table [I] we have collected these resonances with parameters known from Particle Data 
Group book j28|. In the present analysis we shall take into account the following scalar: 
a(600), /o(980), /o(1500), /o(1710), tensor: /2(1270), /^(1525), /2(1565), /2(1950) as well 
as spin-4: /4(2050) resonances. While the tttt decay widths are usually well known, the 77 
decay widths are known only for some of them. Therefore, studying the 77 — >■ tttt data may 
help in extracting the latter quantities. Below we shall present both, theoretical estimates of 
some of them, as well as a fit to experimental total cross sections and angular distributions. 
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TABLE I: Parameters of resonances used in our calculations. 



Resonance 


ruR (MeV) 


Tn (MeV) 


(MeV) 


(keV) 


r r 

i^T-i.. (keV) 


/o(600) [28] 


600 


400 


400 


0.5 




/o(980) [28J 


980 


50 


51.3 [34] 


0.29 




/2(1270) [28] 


1 275 


185 


156.9 


3.035 




/o(1500) [28] 


1 505 


109 


38 




0.033 [29] 


/^(1525) [28] 


1 525 


73 


0.6 


0.081 




/2(i565) m 


1 570 


160 


25 


0.7 




/o(1710) [28] 


1 720 


135 






0.82 


/2(1950) [28j 


1 944 


472 






1.62 


/4(2050) [28] 


2 018 


237 


40.3 


0.7 





In most cases, PDG gives values of the resonance parameters: m^j, F^, F^^ and F 



77- 



We use somewhat smaller values of decay widths for /o(600) resonance than given in the 



PDG book. In the PDG book 28|] a broad range of parameters is given: Fr=(600-1000) 

MeV, F^^ =(1.2-10) keV. In our case we find: Fjj=400 MeV, F^^ =0.5 keV as the best 

fit parameters. In addition, we assume that Br(/o(600) — )■ 7r7r)=100%. For /o(1500) PDG 
r r 

gives T7^7r7r _Q Q33 j^gY ^^ich was obtained by the Belle Collaboration which includes in 
addition in their fit the /o(1370) resonance [1^. Using this value, we obtain: F^^=0.1 keV. 
For /4(2050) the situation is even more complicated. Different values of the ratio (see last 



T^^T 7V7V 



=0.0231 keV, M: 



< 



column in Table H]) has been given in the literature: 



1.1 keV. In our fit to the Belle experimental data (to be shown in the Result section) we 
find: = 0.12 keV. 

i R 

The angular distribution for the s-channel resonances can be written in the standard 
(typical for Feynman-diagrams) form: 



da (77 —7- tttt) 
dz 



E 



4 



— m 



w 
2 



-\M (Ai,A2) 



4 X 647rW2 



(2.11) 



where the helicity-dependent resonant amplitudes must be modelled, as we do not have 
a priori microscopic models of the coupling of two photons to high-spin resonances. Our 
parameterizations for the /o, /2 and /4 resonances are: 



A^(Ai,A2) 



647rW2 X Svr (2 J + 1) (^^) TrTr (W) Br {R 77) Br {R tt+ZOtt-Zo) 



X 



V2S 



Ai,A2 



V|;for /2(1270),/^(1525),/2(1950) 
1^2°; for /2(1565) 
[ Vi';for/4(2050) 



> X exp 



.(2.12) 



The last (form)factor was introduced to correct the resonance form far from the actual 
resonance position where the simple resonance form is incorrect. We expect A^^ to be much 
larger than the resonance width F^j. In practice we will treat it as an extra free parameter 
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to be adjusted to experimental data. This has importance only for broad resonances with 
r > 0.1 GeV. 

The energy- dependent resonance widths are parametrized as: 



R 



9 

—. m± . 

^ =F^ {W, R) 
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(2.13) 



where the function F-^ {W, R) is the spin dependent Blatt-Weisskopf form factor [31 
F^=°(Vr,/2) = 1 
F^=2 {W, R) 



{Rpnf + 3 {Rpnf + 9 



^J=4 



{W,R) 



{Rpy + 3iRpy + 9 

(RPr)^ + 10 iRpnf + 135 (Rpn)^ + 1575 {RpR^ + 11025 
{Rpf + 10 {Rpf + 135 (Rp)^ + 1575 {Rpf + 11025 



(2.14) 



Above pr 



71 r> 9 

-f-mi,p 



— ml and R = 1 fm. 



The decay of /2(1565) into dipions was studied first in Ref. |32| by the Crystal Barrel Col- 
laboration, while its decay into two photons was studied in Ref. H- The /4(2050) resonance 
was not studied so far in 77 — i- nn reactions. 

We consider two simple models of the amplitude for the tensor meson production: 



oc y2,Ai-A2(^,0) • ('5ai-A2,-2 + 5\^-\2,2) (type A), 
« F2,Ai-A2(^,0) • ^Ai-A2,o (type B). 



(2.15) 



We shall call them model A and B, respectively. 

For /2(1270) production the amplitude is dominantly of the type A 27|. We assume the 
same for the /2(1525) and /2(1950) resonances since they are the same nonet partners of 
/2(1270). 

The /2(1565) resonance is not well understood so far. The Belle Collaboration finds in 
this region both components (A, B) in the partial wave analysis H H [sj. We shall try 
both models in order to describe the experimental data, taking the Tt^^ decay width from 



33[ and fitting F^^. 

For the spin-4 resonance again two simple possibilities come into the game: 



■^A^At 



oc Ki,Ai-A2(^,' 



• (5ai-A2,-2 + ^^Ai-A2,2) (type A), 

• '^Ai-A2,o (type B). 



(2.16) 



While the Belle partial wave analysis [22[ 29, l^] suggests the dominance of the type A form, 
in our analysis we find that type B fits better to experimental data. 



77 



vr vr in a simple coupled channel model with p exchange 



Since the cross section for the 77 — )■ tt+tt reaction is much bigger than that for the 
77 — )■ 7r°7r° reaction, even a small coupling between these channels may modify the cross 
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(a) 




FIG. 3: 77 ^ vr+vr- 



section for the 77 — )■ tt^tt*^ reaction. An example of a process which leads to channel coupling 
is shown in Fig|3l 

The contact amplitude (see diagram (a) in Figj3]) can be written as: 

M''i\u\2) = /2eV%(Ai)e.(A2)T:2 + 



p 



X F{sXu)F\K)F\h)F\k2)-^, (2.17) 

(27r) 



and the corresponding t-channel amplitude (see diagram (b) in Figl3]) as: 



4 + A;f — m2 + ie' 



^ ^ 2~~i Otttt^p + P2) ,2 9 I • 

K"^ — + ZL pTUp ^ ^ f^2 ~ ""^^ + ^£ 

X e,(Ai)(/€^ + fcr)e,(A2)(4-A;2^) 

X F(S,t,^)F2(K)F2(A;0F2(A:2)-^. (2.18) 



The form factors that appear in the above formulas are parametrized as: 
• F{s,i,u) = f^p2(^gj \ 



F{x = t, u, ki, k2) = exp 



A4 ' 



F[x = k) = exp 



where A is, in principle, a free parameter. The resulting amplitude strongly depends on the 
value of the parameter. 
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The u-channel amplitude can be obtained from the t-channel amphtude by interchanging 
7T+ and 71- (A^" (Ai, As) = M' (Ai, A2)). 

The corresponding distribution in z = cos 6, when artificially separating the process, can be 
obtained from the amplitudes as usually: 



Ai,A2 



w 
2 



M\Xi, As) + A^"(Ai, A2) + M%Xu A2) 



47r 



4 X 647rW2 



(2.19) 



D. pQCD mechanisms 

In Ref.^ we have concentrated on perturbative mechanisms. In the present paper we 
shall only summarize the main formulae. Again, as shown in FigJH we shall consider simul- 
taneously the n~^Tr~ and tt^tt'^ channels. 




FIG. 4: The Brodsky-Lepage or hand-bag perturbative mechanisms for large-angle 77 — )• tttt 
scattering. 



1. Brodsky-Lepage mechanism 




FIG. 5: Feynman diagrams describing the 77 — t- tttt amplitude in LO pQCD. 



The basic diagrams of the Brodsky and Lepage formalism are shown in FigJSl The 
invariant amplitude for the initial helicities of two photons can be written as: 

M (Ai, A2) = £ dx £ dy 0^ (x, fil) Tt^' [x, y, /i^) 0, (y, /ij) x (t, u) , (2.20) 
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where for fi^ or /i^: ^^^y = min (^x/y,l — x/y^ \J s(l — z^); z = cos9 [TT]. We take the 
hehcity dependent hard scattering amphtudes from Ref . [18| . These scattering amphtudes 
are different for tt+tt^ and n^n'^. The extra form factor in Eq J2.20l was proposed in Ref. (igI] : 



ppQCD (. 
^ req 



U 



exp 



A2 

, reg 



1 — exp 



u — u„ 

A2 

^ reg 



(2.21) 



where tm the maximal kinematically allowed values of t and u. This form factor 

excludes the region of small Mandelstam t and u variables which is clearly of nonperturbative 
nature. In the present analysis we are even more conservative and try to separate the phase 
space into low and high energies and postulate that the pQCD effects show up only above 
a certain energy. In practice we use the following function which smoothly switches off the 
pQCD contribution at low energies: 



FPQCD{s) = 1 - exp 



{s - Aml)^' 

^pQCD 



(2.22) 



Figj6] illustrates the role of the extra form factors discussed above. 




FIG. 6: Comparison of different forms of the extra form factor. We show results for a new 
separation of low and high energy processes given by Eq. (|2.22p and represented by the red solid line 
{ApQCD = 2.5 GeV). In Ref.fsl we used the form given by Eq. (j2.2ip represented in the figure by 
the green dashed line (Aj-eg = 1 GeV). The dotted line shows the BL pQCD cross section without 
any extra form factor. 



The distribution amplitudes are subjected to the ER BL pQC D evolution [35|, l36|. The 



scale dependent quark distribution amplitude of the pion [37 
of Gegenbauer polynomials: 



38| can be expanded in terms 



(p^ (x, /i^) 



A. 

2V3 



6x(l 



X 



1) flr 



n=0 



(2.23) 
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where the expansion coefficients a„ (/i^) depend on the form of the distribution amphtude 
07r Wu and Huang |39j proposed recently a new form of the distribution amphtude: 




x(l-x){l + B X 2 (2x - 1 



\i 8l3^x{l-x) 



-Erf 



\ 8/3% (1 - x) 



(2.24) 



This pion distribution amphtude at the initial scale (/Xg = 1 GeV^) is controlled by the 
parameter B. This simple model better describes recent BABAR data (4o| . 

In FigJT] we compare the results of the Brodsky-Lepage (BL) approach and experimental 
data for different energies. As can be seen from the figures the BL pQCD approach is not 
enough to describe the "high energy" data, especially for the 77 — t- tt^tt^ reaction. 




0.5 0.6 
|cose| 




FIG. 7: Comparison of BL pQCD angular distributions with Belle data |22l. |41| 



2. Hand-bag mechanism 



The hand-bag approach was proposed in Ref.[l9i]. In this approach the only nonzero 
helicity-dependent amplitudes read 

A+- = A-+ = -4naem—R2n (s) . (2.25) 

tu 

The form factors R2n{s) are of nonperturbative nature and are in principle unknown. In 
practice they were fitted to the experimental total (integrated over experimentally measured 
region) cross section for 77 — 7r"'"7r~ [l9|, assuming that the mechanism exhausts the mea- 
sured cross section at high energy. This is not a necessary condition. In our analysis we 
shall relax this rather restrictive assumption. In the present paper we discuss also respective 
angular distributions. We shall use the parametrization proposed in Ref. [ll 



R2. is) = la^ V la, f^y\ (2.26) 



9s V s y 9s 
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FIG. 8: Comparison of hand-bag angular distributions with Belle data 221. |41|. 



The values of a„, n„, and rig are taken from Ref.|19!]. We have averaged the values for 
different energies: a„ = 1.375 GeV^ n„ = 0.4175, a.^ = 0.5025 GeV^ and = 1.195. 

In FiglHlwe show angular distributions calculated in the hand-bag approach together with 
the Belle experimental data 22, 4l||. Only at the highest energies the hand-bag parametriza- 
tion looks reasonable, however, does not describe properly the ratio of the 77 — > tt^vt^ and 
77 — )■ 7r~^7T~ cross section (see Fig lT^ . At somewhat lower energies the hand-bag and ex- 
perimental distributions differ significantly. The onset of the hand-bag mechanism was not 
discussed so far in the literature. We shall return to the problem in the Result section where 
other mechanisms will be discussed too. 



E. Equivalent photon approximation for AA — t- AAtttt 

How to calculate the production of pairs of particles in nucleus-nucleus collisions was 
explained in detail in Ref. [sf. Here we only sketch the formalism. The total nuclear cross 
section can be expressed by folding the 77 — )■ vrvr subprocess cross section with equivalent 
photon fluxes as: 

a{AA AAtttt; saa) = ^ <3-(77 -> vrvr; Vr^^)^2^,(b) 

X N{ui,hi)N{u2,h2)d%id%2dujidu2 . (2.27) 

In principle, the absorption factor S'abs(b) can be calculated in different models. At not too 
high energies the Glauber approach is a reasonable approach. Here we take a somewhat 
simpler approach and approximate the absorption factor as: 

5L(b) = ^^(b - 2Ra) = e{\h^ - h2\ - 2Ra) . (2.28) 

This form excludes the situation at high energies when the colliding nuclei overlap which at 
high energies unavoidably leads to their breakup. 

In order to simplify our calculation we shall use the following transformations: 

= ^e±^-% duj^du2 = ^dW^^dY^^ . (2.29) 
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Then formula f l2.27p can be rewritten as: 

a{AA AAtttt; saa) = J ^{ll ^ vrvr; Wyy)Slf„{h) 

X N{LOi,hi)N{uj2,h2)^^d'^hid%2dW^^dY^^ . (2.30) 

While in the above approach one can easily calculate only the total nuclear cross section, 
distributions in rapidity of the pair of pions, and the invariant mass of the dipions (see e.g. 
[1, 111), experimental constraints can not be easily imposed. 

If one wants to calculate in addition kinematical distributions of each of the individual 
particles (transverse momentum, rapidity, pseudorapidity) , or impose corresponding exper- 
imental cuts on kinematical quantity of each of dipions, a more complicated calculation is 
required. Then instead of the one-dimensional array (7(77 — )■ 7r7r;VFj), a two-dimensional 
array '^"'^'^^^"""^^ (Wi, Zj) has to be calculated, stored and passed to the code calculating cross 
sections for nuclear collisions. Then an extra integration m z = cos 6 is required in Eqs. fl2.27p 
or fl2.30p . which makes the calculation much more time-consuming. 

Next four-momenta of pions in the vrvr center of mass frame are calculated: 

= Vs/2. , 



^ 2 



Pi = zpn . (2.31) 
The rapidity of each of the pions can be calculated easily as: 

Vi — Y^TTTT ~\~ Vi/mT ; (2.32) 

where yi/T^n = yi/TrniW, z) is the rapidity of one of the pions in the recoil vrvr system of refer- 
ence. The transverse momenta of pions in both frames of reference, as Lorentz invariants, 
are the same. Other kinematical variables are calculated by adding relativistically velocities 

m 

(2.33) 

where ^ 

VZ = ^ (2.34) 

and from the energy-momentum conservation: 

= uj^-uj2, (2.35) 
the energies of photons can be expressed in terms of our integration variables as: 

uji = -^exp(y) , 

CU2 = ^exp(-y). (2.36) 
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Now the pion velocities can be converted to four-momenta of pions in the overall nucleus- 
nucleus center of mass frame. Then angles or pseudorapidities of pions can be easily calcu- 
lated and experimental cuts can be imposed in addition. 

Distributions in kinematical variables are obtained by corresponding binning into his- 
tograms. A smooth distributions in transverse momenta or pseudorapidities require many 
points in W and z. If only experimental cuts are imposed the number of points can be 
smaller. In the present approach we are interested first of all in two-pion invariant mass 
distribution but including cuts of real experiments. 



F. AA AAp 



In nucleus-nucleus collisions another mechanism exists which contributes to the vr+vr" 
channel - coherent production of p° mesons and its subsequent decay into the tt^tt^ channel 

quantum Glauber , Glauber-Gribov |45|, 146 



This mechanism was studied theoretically in the classical Glauber 

and Color Glass Condensate 14 



approaches. 



The rapidity distribution of can be written in an almost model-independent way as: 



da dNr^iuji) dN^{u2) 

-r = 1^1 J <^fA^VA + <^2 1 O-A-y- 

dy duji duj2 



yAV 



(2.37) 



Essentially the models enter in calculating the cross section for 7^4 — )■ VA processes. While 



for heavy quarkonium production the situation is more transparent [481], for light meson 
production the situation is somewhat less certain. We shall not discuss the issue here, we 
rather refer the reader to the existing literature (see e.g. 49|). In the most rudimentary 
approach: 



a 



tot\ 
pAj 



dtFlit) 



(2.38) 



where cr*'^^ 



T^R\ and tr. 



III. RESULTS 

Before we go to the nuclear processes we wish to discuss individual subprocesses 77 — )■ 
7r+7r~ and 77 — )■ -k^-k^ which are interesting by themselves. 



77 



vrvr 



We wish to show first cross section for 77 — )■ tt+tt^ and 77 — )■ vr'^vr'^ processes as a function 
of the subprocess energy. We shall include the pion exchange mechanism for the 77 — )■ tt+tt" 
reaction, several resonance states, pion-pion rescattering important, as will be shown below, 
for 77 — 7- Ti^Ti^ reaction where in the leading order the pion exchange continuum vanishes. 
At somewhat higher energies we include also the Brodsky-Lepage pQCD and hand-bag 
mechanisms. In FigJH] we show our model calculations against world data for 77 — )■ tttt 



22I . l34i l41l |50H56[ . We get a good agreement with all available data points for the first time 
in such a large range of energies. This makes our hybrid model for 77 — t- tttt well suited for 
the predictions of the cross sections for nucleus-nucleus collisions. 
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FIG. 9: Results of the fit: the total cross section (77 — t- tt'^-k for | cos9\ < 0.6, 77 — )■ vr'^Tr*' for 
I cos 6*1 < 0.8). 




FIG. 10: Contributions of different s-channel resonances to the 77 — t- tt+tt reaction. 

Our model amplitude is a sum of many ingredients as discussed in the previous section. 
In FigHniand FiglTT]we show some somewhat artificially ^ separated ingredients to the final 
spectrum. We show s-channel resonance (Fig JTOj) and continuum (Fig JTT]) contributions sep- 
arately. At low energies, the cr(600) s-channel resonance as well as pion rescattering due to 
p-meson exchange play crucial role. At somewhat higher energies the spectra are dominated 
by the contribution of the tensor /2(1270) meson. At still higher energies the situation is 
less clear. The smallness of the cross section at W ~ 1.5 GeV can be understood as an 



We add up amplitudes not cross sections. 
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W,, = M,,[GeV] 



W,, = M,JGeV] 



FIG. 11: Contributions of different continuum processes. The meaning of the curves is explained 
in the figures. 

interference of right wing of the tensor /2(1270) resonance and other resonances (/2(1525), 
/o(1500), /2(1565)). At around W ^ 2 GeV, we observe a new irregular behaviour. This 
can be understood as the presence of s-channel spin-4 contribution. In order to describe 
the angular distribution we use a simple model assuming that the spin-4 amplitude is pro- 
portional to A^AiA2 ~ ^,Ai-A2^Ai-^ (type B). This is similar to what was found from a 
phenomenological analysis in Ref.[22|. A much better description of the data is obtained 
when the /4(2050) resonance is included. The corresponding radiative decay width is treated 
as a free parameter. We get Tf^^^^ = 0.7 keV ± 0.1 keV. This is rather a rough estimate 
than a rigorous statistical analysis. The problem is in the background which is in this range 
of invariant masses rather difficult to control. Our analysis is a first trial to extract the 
radiative decay width of the /4(2050) resonance. The found value is much smaller than 
that for the /2(1270) (3.035 keV) and comparable to that for /o(980) (0.29 keV) (see Table 
H]). We shall not discuss consequences of better understanding of the structure of /4(2050) 
resonance. This certainly goes beyond the scope of the present paper. 

Clearly, the tensor /2(1270) is the dominant resonance. Other resonances give much 
smaller but sizeable absolute cross section. However, their separation is not an easy task 
especially from the analysis of total cross section. As was discussed recently [13], the cross 
section for the production of glueballs (G) in the nucleus-nucleus AA — )■ AG A reaction could 
be large. However, selecting a final state is in this context a very important issue. The tttt 
final state is an obvious candidate. Our analysis shows that searching for glueballs in the vrvr 
channels may be extremely difficult, if not impossible. This makes the heavy ion processes 
not the best tool to study glueballs or other exotic mesons. We think that a precise fit of 
the present Belle angular distribution is probably better than any fit to the future nuclear 
data. 

In Figim we show the contributions of the Brodsky-Lepage and hand-bag mechanisms. 
Please note that our normalization of the hand-bag mechanism was adjusted to describe 
high-energy experimental data together with the Brodsky-Lepage mechanism. In Ref. lol 



only the hand-bag mechanism, ignoring the Brodsky-Lepage mechanism, was fitted to the 
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high-energy data. For the vr"*"??" channel we show in addition the continuum distribution 
corresponding to pion exchange for two different values of slope parameter and for the 
tt^tt" channel rescattering contribution due to t-channel exchanges. The corresponding 
contribution falls quickly with increasing energy and plays an important role only up to 
W = 1 GeV. In FiglH] (right panel) we show the sum of resonances, BL pQCD and hand-bag 
mechanisms - lower line. The upper curve includes in addition the exchange contribution. 
Adding the p exchange contribution allows us to better describe experimental data for small 
invariant masses. 




|cose|; W=0. 8025 GeV |cose|; W=1. 1525 GeV |cos9|; W=1 .4975 GeV 





FIG. 12: Examples of angular distributions for the 77 — )• tt+tt reaction. 



In FigHJ] and [13] we compare our predictions to the angular experimental data measured 



by the Belle Collaboration. Experimental data points are taken from Refs. [34, |41| for 
77 — )■ 7r+7r~ and from Ref. j22| for 77 — )■ 7r°7r° reactions. In Fig. [12] (the first four panels) 
we show the sum of contributions of many mechanisms for two different values of the slope 
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FIG. 13: Examples of angular distributions for the 77 — >• tt ir reaction. 



parameter in the hadronic form factors (see below Eq. (12. Sp ). We show how the final result 
depends on the value of this parameter. In Ref. the results for i? < 4 GeV were shown, 
but here we obtain a better description of the data when we use the exponential form factor 
with B^Tj=6 GeV~^. The description for W =1.4975 GeV shows that one should include also 
D wave, but we have no theoretical guidance how to do it. For larger W, the contribution 
from pQCD mechanisms starts to dominate. In Fig {13] we show the angular distribution for 
77 — > 7r°7r° process for nine different energies ranging from 0.61 GeV to 4.05 GeV. Here 
we use the same parameters of resonances (masses, widths, phases) as for the 77 — )■ vr+vr" 
reaction. Here we have a small problem with good description for two energies (1.05 and 
2.34 GeV). In addition, we show the cross section corresponding to the sum of processes 
(amplitudes) with and without contribution of the p^-meson exchange. We can observe 
that the inclusion of this contribution is important only for W < 0.6 GeV (see. FiglHl right 
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FIG. 14: The ratio of the 77 — )• vr^vr^ and 77 — t- vr+vr cross sections as a function of W- 



77 



panel). In conclusion, we have obtained a reasonable description of the angular experimental 
data, which results in excellent description of the total cross section (see Figj9]). 

In Fig{T3]we present the ratio of the neutral and charged pion pair cross sections. The 
upper line represents the hand-bag model result, which is independent of 6 and is equal 
to ^. Our result, which includes the BL pQCD and hand-bag contribution simultaneously, 
describes the experimental data measured by the Belle Collaboration 41|. We have a 
problem with correct description of the data at W = (2.4-2.7) GeV. We think that this is 
due to not exact description of angular experimental data for 77 — )■ tt^tt" reaction in this 
range of energy. However, our result shows that for the correct description of high energy 
cross section it is necessary to take into account also the BL pQCD contribution. 

In FigHSl we summarize our elementary cross sections as a function of ^/s^^ and cos6'. 
Here transverse momentum cuts were imposed as explained in the figure caption. The lower 
cut roughly corresponds to the lower cut of the ALICE experiment. It is not clear for 
us if such a cut is sufficient to regularize the calculations of the BL pQCD and hand-bag 
mechanisms. Therefore we present also a result with slightly bigger cut. Another option 
would be to cut in cos 6* as is usually done at e+e~ colliders. The irregularities observed at z 

± 1 are caused by small number of bins in the figure. Such two-dimensional distributions 
as shown in Fig|T5] are then used in nuclear calculations discussed in the next subsection. 

B. PbPb PbPhTT-K 



In order to check our nuclear calculation, in Fig{T6] we show interesting distribution 
in impact parameter between the two lead nuclei. The distribution is different from zero 
starting from the distance of 6 = 2Ra ~ 14 fm and extends till " infinity" . This means that 
a big part of the cross section comes from the situations when the two ^^^Ph nuclei fiy by 
very far one from each other. This figure clearly demonstrates the ultraperipheral character 
of the discussed process of exclusive dipion production via photon-photon fusion. 

The total cross sections for ^Jsnn = 3.5 TeV are collected for both vr+vr" and -k^-k^ 
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FIG. 15: -^{W,z) for 77 — t- vr^vr (left) and 77 — t- vr'^vr'^ (right) for two different cuts on pion 
transverse momenta: > 0.2 GeV (top) and pt,n > 0.5 GeV (bottom). 



TABLE II: Cross sections for different lower cuts on pion transverse momenta at ^/sJfN = 3.5 TeV. 



Pt,min [GeV] 



TT ' TT 



[mb] 



mb] 



0.2 
0.5 
1.0 



46.7 


8.7 


12.1 


5.1 


0.08 


0.05 



channels for different lower cuts on pion transverse momentum in Table [TTl 

The distribution in invariant dipion mass shown in Fig{T7] is the most interesting one. 
We show somewhat theoretical distribution for the full phase space (upper solid line). We 
show distributions for the PbPb — )■ PbPbir^n' (left panel) and for the PbPb — )■ PbPbT^^T^^ 
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FIG. 16: The impact parameter distribution for the PbPb — t- PbPb-K^ir and PbPb — t- PbPbTT^n* 
reactions at ^snn = 3.5 TeV. 




0123456 01234 



M,., [GeV] M„„^„[GeV] 

FIG. 17: Total cross section for 77 — J- t\^-!\~ (left panel) and 77 — )• Tr^vr" (right panel) as a 
function of the system energy at the LHC energy y^sjviv = 3.5 TeV and -pt.-K > 0.2 GeV. We show 
in addition results with extra cuts on z. 

(right panel) reactions for the full phase space (upper lines), for | cos6'| < 0.9 (middle lines) 
and for the angular range corresponding to the experimental limitations usually used for 
the 7171 production in e~^e~ collisions (| cos6'| < 0.8). At lower energies {W < 1.5 GeV) the 
result does not depend on the angular cuts. The big differences start in the region where 
the elementary cross section is described by the BL pQCD and hand-bag mechanisms. 

FigJTSl shows two-dimensional distributions in pseudorapidity of charged pion (left panel) 
or pseudorapidity of neutral pion (right panel) and transverse momentum of one of the pions. 
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FIG. 18: -3;^^$^ for the ^"«P6 + ^ 
208p5 + 208 p5 4" 208p5208p^^o^o reaction (rij 

and pt,^ > 0.2 GeV. 



208p5208p5^+^- reaction (left panel) and the 
;ht panel) at the LHC energy ^/snn = 3.5 TeV 



With higher pt^^ values, the pseudorapidity distribution becomes somewhat narrower. 




FIG. 19: ^ and |^ (left panel) and ^ (right panel) for the LHC energy ^/i^ = 3.5 TeV. 



Let us start now presentation of theoretical differential distributions which can be in 
principle measured. In the left panel of FiglT9]we show distributions in rapidity of individual 
pions from the pair (solid line) as well as distribution in pseudorapidity of the same pion 
(dashed line). Right panel illustrates nuclear cross section as the function of the pion pair 
rapidities. Here we have imposed extra cuts on pion transverse momenta: pt,-K > 0.2 GeV 
(solid lines) and pt^n > 0.5 GeV (dashed lines). In addition, we compare the nuclear cross 
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section for vr+vr (upper curves) and for vr^vr'^ (lower curves) production. These distributions 
(in tItt and F,r7r) look fairly similar. 




0.5 1 1.5 2 2.5 



FIG. 20: ^ of one of pions for the LHC energy = 3.5 TeV and pt,^ > 0.2 GeV. 

The distribution in transverse momentum of one of the pions is shown in Fig. [20] for the 
full phase space for PbPb — )■ PbPbir'^Tc" (solid line) and for PbPb — )■ PbPbTT^ir^ (dashed 
line), respectively. The distribution shown in Figj20]is for the full range of (pseudo)rapidities. 
The ALICE detector can cover only a part of the whole rapidity range. 




FIG. 21: "'^^ _ for the 208p6 + 208p^ ^ 208 ^^208^^^+^- reaction at the LHC energy ^/I^ = 
3.5 TeV and pt,^ > 0.2 GeV. 

In Figl2T]we show a two-dimensional distribution of pseudorapidities of pions. The cross 
section is concentrated along the diagonal rji = 772- The range covered by the ALICE 
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detector (a square (-0.9,0.9) x (-0.9,0.9) centered around r/i = 772 = 0) is rather small. The 
pseudorapidity distribution of one of the pions is shown in FigJT9] (left panel) by the dashed 
line. 



C. p^{770) and p'^(1450) photoproduction in PbPb collisions 



The TT+TT pairs produced in ultraperipheral photon-photon (sub) collisions constitutes 
a background for, e.g., exclusive production of the p° meson. This reaction was studied 



experimentally at RHIC [57[. Theoretical calculations for LHC energies predict very large 



cross section of the order of a few barns (see e.g. [43|, |44, |46|, |47[), i.e., more than two 
orders of magnitude larger than the contribution of the mechanism discussed here (see also 
Table UTTll . The p° spectrum is, however, concentrated in pion-pion invariant mass around 
the p-resonance position. It is an open question how far from the resonance position one 
can trust a rather naive Breit-Wigner form. In photoproduction of meson on protons one 
often fits the following form: 



da{PbPb PbPbp^ 
dMZ 



A„o — — 
P M2 



M^7r^pOr(M7r 



+ A 



(3.1) 



where ApO is the amplitude for the Breit-Wigner function, A^^r is the amplitude for the direct 
TT+TT" production and T{Mt^t^) is the momentum- dependent width of the p° resonance. In 



our presentation below the values of parameters are taken from Ref. [5f 
total cross section was assumed to be 5 barn (see Table UTTll . 



and the integrated 



TABLE III: Cross sections for coherent p and p' production for LHC. 



Reference / model 


^SNN [TeV] 


c^CF''' [b] 


-i:/""^ [b] 


KN [43] 


5.6 


5.2 




FSZ [4^ 


5.5 


9.538 


2.216 


RSZ [46] 


5.5 


9.706 




IKS [62] /KST-R 


5.5 


4.9 




/KST 




4.36 




/GBW 




3.99 




/VDM 




10.03 




CM [63] 


5.5 


10.069 





In Figj22]we show the contribution of the p°(770) meson integrated over the whole phase 
space, without any extra cuts. While the predictions for the production are rather uncer- 
tain (within factor two or three), our predictions for the 77 process should be fairly precise 
(a few percent precision). A comparison with real data should be therefore interesting. At 
larger dipion invariant masses, other radial excitation may appear in addition. The p''(1450) 
state is a first candidate of this type. Since the continuum single-photon photoproduction 
7r"'"7r~ background in this region was not studied, we have no information about its inter- 
ference with the p''(1450) resonance. Therefore in the following we shall show only the 
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resonance contribution and leave the question of the background for future analyses, includ- 
ing perhaps the ALICE experimental data. Here we parametrize the p'^(1450) contribution 



as: 



da{PbPb PbPbp^{U50)) 



^p0(1450)^ 



^M^^mpO(i45o)r(il4 



r(M. 



(3.2) 



The cross section for the photoproduction of p'^(1450) was calculated only in Ref. 4^]. This 
calculation finds that the cross section for the p°(1450) resonance is about a factor 5 smaller 
than for p°(770). Here we are discussing the signal in the tt+tt^ channel, therefore a corre- 
sponding branching fraction has to be included in addition. Unfortunately this branching 
fraction is not well known. The Crystal Ball Collaboration has measured only the ratio of 
the cross section for two and four pion channels [59i]. They have found i?r(27r)/i?r(47r) = 
0.37 ± 0.1. In this situation we can only calculate an upper limit for the two-pion branching 
fraction: 



Br (p°(1450) TT+TT" 



P 



2-7T 



P2-K + Pi-K + Pc 



< 



P 



2lT 



P2JP 



An 



other 



P2n + PiiT P2iv/ Patv + 1 



0.27 . (3.3) 
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FIG. 22: Invariant mass distribution of vr^vr from the decay of /)°(770) and /j°(1450) photo- 
production (resonance contributions represented by the dashed and dotted lines) and our 77 fusion 
(soUd Une) in ultraperipheral Pb-Pb collisions at y/s^N = 3.5 TeV. 

In Figl22] the resonance contribution for the p°(1450) is 5/0.27 times smaller than for 



p°(770) as suggested by the calculation 4^ and the above upper limit for the branching 
fraction into pions. It is therefore clear that this will be an upper estimate for the p'^(1450) 
contribution. The relative contribution of p°(1450) is still somewhat larger than observed 
for instance in the electroproduction of these states at HERA [60]. In Ref. [6.l| an anomalous 
(strongly nucleus mass and photon virtuality dependent) production of p°(1450) has been 
predicted for incoherent production in the 7*^4 — )■ p°(1450)X reaction. We wish to emphasize 
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in addition that our estimate does not take into account experimental cuts of any concrete 
experiment. 

The photoproduction contribution cross sections are above our photon-photon fusion one. 
It is not clear in the moment how the kinematical cuts may change the proportions of the 
two mechanisms. This will be studied elsewhere. 

In contrast, the PbPb — )■ PbPbTi'^Tr^ reaction is free of the photoproduction mechanisms. 
However, it is not clear for us if, in this measurement is possible. 



IV. CONCLUSIONS 

In the present paper we have discussed first 77 — >■ 7r°7r° and 77 — >■ tt+tt" reactions and we 
have shown that different mechanisms contribute to such processes below ^ys^ < 6 GeV. 

These are: pion exchange for the tt+tt" channel, several resonances such as /o(600), 
/o(980), /2(1270), /o(1500), /^(1525), /2(1565), /o(1710), /2(1950) and /4(2050) and pQCD 
mechanisms proposed first by Brodsky and Lepage as well as hand-bag mechanism proposed 
by Diehl, Kroll and Vogt. We have included also pion-pion rescattcring which leads to a 
coupling between the tt+tt" and vr'^vr*' channels. We have shown that the inclusion of such 
processes can help in understanding energy dependence of the 77 — )■ nn processes as well as 
angular distributions only at very low energies. We have not tried to get a perfect fit to the 
data as several details of the model amplitudes are not known. We have slightly adjusted 
some parameters to get a reasonable description of the data. 

We have found that the decay width (7(600) — )■ 77 is much smaller than that found in 
other partial wave analyses. The data (total cross section as well as angular distributions) 
close to the position of /2(1270) peak can be described assuming the dominance of the 
amplitude 

OC Y,,,,_,,{e, <t>) ■ (5m-A„-2 + <^Ai-A„2) . 

We have included also the tensor /2(1565) resonance, considering two models of the corre- 
sponding 77 — > TTTT amplitude. A better description of the data has been obtained using 

i.e., space, different than for /2(1270). This may mean that the structure of this resonance 
is completely different than that for /2(1270), /a (1525) or /2(1950) {qq states). 

Inclusion of the /4(2050) spin-4 resonance improves the spectrum and angular distribu- 
tions around the resonance position. By fitting to the Belle Collaboration data we have 
found (2050) ->77 = 0.7 ± 0.1 keV. 

At still higher energies the situation is even less clear. While in the tt+tt" channel the 
Brodsky-Lepage mechanism with distribution amplitude discussed recently in the context 
of the BABAR data describes the data, in the "smaller" ti^ti^ channel more processes may 
contributed in addition to the BL pQCD mechanism. Another possibility is the hand-bag 
mechanism proposed some time ago by Diehl, Kroll and Vogt. The hand-bag mechanism 
may be important only at energies W > 3 GeV. 

At high energies the angular distributions show an enhancement at large | cos^| for the 
77 — 7- vr'^Tr^. Such an enhancement is predicted both, by the Brodsky-Lepage and by the 
hand-bag approaches. A proper mixture of both processes provides a better description of 
the experimental data than each of them separately, as done usually in the literature. 
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Having fixed the details of tlie amplitude we have calculated total cross sections and 
angular distribution as a function of 77 energy for both 77 — )■ tt+tt" and 77 7r°7r° 
processes. These two subprocess energy- dependent cross sections have been used next in 
Equivalent Photon Approximation in the impact parameter space to calculate for the first 
time corresponding production rate in ultraperipheral ultrarelativistic heavy ion reactions. 
In this calculation we have taken into account realistic charge distributions in colliding 
nuclei. 

We have calculated both total cross sections at LHC energy, as well as distributions in 
rapidity and transverse momentum of pions and dipion invariant mass. The calculation 
of distributions of individual pions is slightly more complicated in the b-space EPA. The 
distributions in dipion invariant mass have been compared with the contribution of exclu- 
sive p° — )■ TT^TT^ production in photon-pomeron (pomeron-photon) mechanism taken from 
the literature. Close to the p° resonance the 77 — )■ tt+tt" mechanism yields only a small 
contribution. The 77 contribution could be, perhaps, measured outside of the p° resonance 
window. A detailed comparison with the absolutely normalized ALICE experimental data 
should allow a quantitative test of our predictions. Imposing several experimental cuts may 
enhance the 77 — )■ vrvr contribution. 

In the present paper we have analysed exclusive production of two pions. In principle 
exclusive production of four pions could be also interesting. First experimental study has 
been presented in js^]. Interesting theoretical study has been performed in (63 |. 



Appendix A: A comment on partial wave decomposition 

In our model, the amplitude for the 77 — )■ vrvr process is a coherent sum of different 
contributions (amplitudes) for different mechanisms: 

-Ma„a.W=EM'!a.W, (A1) 

k 

where Ai, A2 are initial photon helicities. There are only four combinations of helicities 
defining the, in general, complex amplitude: (-1,1), (1,-1), (1,1), (-1,-1). 

As a consequence the differential distributions for the 77 — )■ vrvr reaction can be written 
in terms of spherical harmonics numbered with pion-pion angular momentum projection 
(M = A1-A2) 

^ = |Co,o(Vr)Fo,o(^,0) + C2,o(W^)>^2,o(^,</') + C4,o(TV)n,o(^,0) + (-)l' 
ail 

+ 1^^2,2(1^)1^2,2(^,0) + C4,2(i^)n,2(^,0) + (..or 

+ \C2,-2mY2,-2{z,(p) + C,,.2iW)Y,,.2iz,(P) + {...)\' . (A2) 

The expansion coefficients Cl,m are strongly energy dependent. In general, resonance con- 
tributions (L = 0, 2, 4) interfere with slowly-energy-dependent backgrounds (L = 0, 2, 4,...) 
of different origin discussed in previous sections. 

Recently the Belle Collaboration (22I . 29 1 has performed a simplified fit to the experimental 
angular distributions of the type: 

^ = |Co,o(Vr)ro,o(z)|' + |C2,o(l^)r2,o(^)|'+|C4,ol4,o(^)|' 

+ 1^2,2(1^)^-2,2(^)1' + |C4,2(W^)n,2(^)|' . (A3) 
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Our model shows that such a simphfied partial wave expansion has not always model jus- 
tification and may therefore lead sometimes to unphysical solutions. However, the correct 
partial wave decomposition (see (Eq JA2p ) of experimental angular distributions is in practice 
too difficult, if not impossible. 
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